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Abstract 

Canonical  forms  and  canonical  variables  for  inviscid  flow  problems  are  derived.  In  these  forms 
the  components  of  the  system  governed  by  different  types  of  operators  (elliptic  and  hyjierbolic) 
are  separated.  Both  the  incompressible  and  compressible  cases  are  analyzed  and  their  similarities 
and  differences  are  discussed.  The  canonical  forms  obtained  are  block  uppe-  '  riangular  operator 
form  in  which  the  elliptic  and  non-elliptic  parts  reside  in  different  blocks.  The  full  nonlinear 
equations  are  treated  without  using  any  linearization  process.  This  form  enables  a  better  analysis 
of  the  equations  as  well  as  better  numerical  treatment.  These  forms  are  the  analog  of  the 
decomposition  of  the  one  dimensional  Euler  equations  into  characteristic  directions  and  Riemann 
invariants. 
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1  Introduction 


In  the  past  decade  a  substantia]  effort  has  been  invested  in  numerical  solutions  of  the  Euler  equations 
in  one  dimension.  These  were  based  on  characteristic  decomposition  of  the  flow,  i.e,.  a  decomposi¬ 
tion  into  three  subsystems  each  of  which  correspond  to  a  simple  propagation  along  an  appropriate 
line  in  space-time.  Efficient  discretization  schemes  were  based  on  this  representation.  We  refer  to 
[2]  for  extensive  reference  on  the  subject.  This  approach  is  usually  extended  to  multidimensional 
problems  by  splitting  it  into  a  sequence  of  one-diinensional  problems. 

The  above  approach  is  limited  and  does  not  capture  the  real  structure  of  the  solution  for 
multidimensional  flow  problems.  Moreover,  the  analysis  is  usually  based  on  the  time-dependent 
problem  although  the  case  of  interest  is  the  steady-state  problem.  Consequently,  an  important 
feature  of  the  system  such  as  its  mixed  nature  has  not  been  given  proper  attention.  This  paper  is 
concerned  with  multidimensional  steady-state  inviscid  flow  problems. 

An  analysis  which  is  analogous  to  the  decomposition  of  the  one  dimensional  Euler  equations 
into  characteristic  directions  and  the  Riemann  invariants,  is  presented  for  multidimensional  steady- 
state  inviscid  flow  problems.  Thi*^  methodology  can  be  applied  to  more  genera]  partial  differential 
equations,  such  as  the  steady  state  and  time  dependent  viscous  flow  equations.  This  analysis  of 
the  equations  break  up  the  full  problem  into  its  fundamental  (irreducible)  subprcblems  and  also 
describe  the  interaction  between  them.  It  leads  to  a  better  understanding  of  the  equations  and  to 
new  type  of  schemes  that  better  represent  the  physical  behavior. 

The  analysis  is  presented  for  incompressible  and  compressible  Euler  equations  both  both  in 
two  and  three  dimensions.  It  shows  in  a  clear  way  the  existence  of  quantities  that  propagate  along 
streamlines,  and  other  quantities  whose  behavior  is  directionally  unbiased  for  subsonic  compressible 
flows  or  incompressible  flows.  Supersonic  flows,  on  the  other  hand,  both  in  two  and  three  dimensions 
are  hyperbolic  with  respect  to  the  stream  direction.  A  change  of  variables  is  introduced  in  order 
to  bring  the  system  into  its  canonical  form,  and  the  resulting  new  variables  are  called  canonical 
variables.  The  canonical  forms  are  block  upper  triangular  form,  where  the  diagonal  blocks  consists 
of  the  basic  components  of  the  systems  and  off  diagonal  ones  represent  some  interaction  between 
the  different  subsystems. 

For  two  dimensional  incompressible  Euler  equations  the  canonical  variables  are  velocities  and 
total  pressure  p  -f  q^/'2  and  in  three  dimensions  they  are  velocities,  total  pressure  and  normalized 
helicity,  namely,  (a>  ■  V)/(V  •  V).  Canonical  variables  other  than  the  velocities  propagate  along 
streamlines  and  form  the  non-elliptic  part  of  the  system.  The  velocity  components  satisfy  an  elliptic 
system  with  a  forcing  depending  on  the  hyperbolic  parts. 

For  the  the  two  dimensional  compressible  Euler  equations  the  canonical  variables  are  velocities 
and  total  enthalpy,  and  in  three  dimensions  a  normalized  helicity,  namely,  (u?  •  pV)/(pV  •  pV)  is 
the  additional  canonical  variable.  Quantities  other  than  the  velocities  propagate  along  streamlines. 
The  velocities  are  governed  by  an  elliptic  system  for  subsonic  regimes  and  by  a  hyperbolic  system  in 
supersonic  regimes.  The  structure  of  the  canonical  forms  shows  in  a  clear  way  the  relation  between 
the  full  problems  (incompressible  and  compressible  Euler  equations)  and  simplified  models  that 
have  been  studied  in  the  past,  such  as  (div,curl)  systems  in  two  and  three  dimensions,  and  their 
compressible  versions.  Canonical  boundary  conditions  are  discussed  in  this  framework.  These  are 
the  analog  of  characteristic  boundary  conditions  for  one  dimensional  compressible  flows. 


1 


2  On  Canonical  Forms  for  Systems  of  Partial  Differential  Equa¬ 
tions 

Differential  equations  of  constant  coefficient  are  reduced  by  the  Fourier  transform  to  algebraic 
equations  in  terms  of  symbols  which  are  in  general  matrix  polynomials  in  several  variables.  Let 
P{D)  be  a  differential  operator  of  the  form 


P{D)=Yla^D‘^  (2.1) 

|a|<m 

where  D  =  (D^,. . Dj),  a  =  (a\ , . . . ,0^)  and  Z?"  =  ’  . . . Dj"*  is  the  usual  multi-index  notation. 

The  Ua  are  I  x  I  real  valued  matrices.  We  want  to  study  solutions  to  the  problem 


P(D)u(x)  =  f(x)  X  eil  (2.2) 

B(D)u{x)  =  g{x)  X  e  d^l  (2.3) 

where  D  is  a  domain  in  and  B(D)  is  an  appropriate  boundary  condition  to  be  specified  later. 

Fourier  analysis  in  full  and  half  space  are  carried  out  in  order  to  analyze  the  interior  and  bound¬ 
ary  properties  of  the  solutions,  respectively.  The  symbol  associated  with  P{D)  is  the  polynomial 


P(x,tO=  S 

|a|<m 


(2.4) 


where  (  =  (^1 , . . . , ^<j)  and  and  its  use  transforms  the  constant  coefficient  differential 

equation  into  the  algebraic  equation 


HiOm  =  Ho-  (2.5) 

For  elliptic  systems  where  the  matrix  polynomial  satisfies 

|detP(zOI>C|e|""*  (2.6) 

(2.7) 

extensive  research  has  been  done  with  regard  to  smoothness  properties  of  the  solution  including  the 
effect  of  boundary  conditions.  The  case  of  hyperbolic  systems  of  differential  equations,  i.e.,  when 
the  principal  part  of  detP(i^)  has  real  roots  for  ^  9^  0,  has  been  studied  in  the  context  of  time 
dependent  problems. 

For  the  appbcations  discussed  in  this  paper  one  is  interested  in  analyzing  steady-state  problems 
which  are  mixed  elliptic-hyperbolic.  That  is,  part  of  the  system  is  governed  by  hyperbolic  equations 
(with  respect  to  some  direction  in  space)  and  the  rest  of  the  system  by  an  elliptic  sub.system,  and 
some  interaction  exists  between  the  two  subsystems. 

The  first  step  in  the  analysis  of  a  mixed  system  is  the  identification  of  the  different  subsvstems, 
i.e.,  elliptic  versus  hyperbolic.  Next  comes  the  identification  of  the  variables  corresponding  to 
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different  type  of  behaviors,  followed  by  a  reduction  of  the  system  to  block  triangular  form.  The 
diagonal  blocks  are  the  ba.sic  building  blocks  of  the  system. 

We  discuss  two  types  of  reduction  to  canonical  forms.  In  the  first  case  the  new  set  of  variables 
are  linear  combination  of  the  old  ones,  via  constant  matrices.  In  the  other  case  a  linear  combination 
via  polynomials  is  used.  Thus,  the  new  variables  in  the  second  method  involve  linear  combinations 
of  the  old  ones  and  their  derivatives. 

C-Reducibility 

Assume  that  there  e>iSt  invertible  matrices  S,T,  (independent  of  such  that  SP{^}T  is  of 
block  upper  triangular  form,  i.e.. 


SP{^)T  = 


/  PniO 


PuiO 

P22{^) 


PkkiO  / 


=  A 


(2.8) 


where  the  subsystems  PvJ  =  1 , . . . ,  A:  are  irreducible. 

Taking  the  determinant  of  both  sides,  and  assuming  that  det(5)det(T)  =  1,  we  obtain 


det  PiO  =  det  Ai(0  det  P22{0  •  •  -det  Pu(0  (2.9) 

The  blocks  on  the  diagonal,  correspond  therefore,  to  the  primitive  building  blocks  of  the  system, 
and  the  off  diagonal  ones  represent  the  interaction  of  the  different  basic  blocks. 

Definition  1  A  matrix  polynomial  P{^)  will  be  called  C-reducible  if  there  exist  invertible  matrices 
S,T  independent  of  ^  such  that  A(^)  =  SP{^)T  is  of  block  triangular  form. 

We  will  refer  to  this  form  of  the  matrix  polynomial  as  the  canonical  form.  It  induces  a  canonical 
form  for  the  differential  operator  via  the  Fourier  transform. 

Thus,  we  have  the  following  theorem. 

Theorem;  A  necessary  condition  for  an  I  xl  matrix  polynomial  P(i^)  to  be  C-reducible  to  a  block 
triangular  form  is  that  its  determinant  admits  a  factorization  into  lower  order  polynomials. 

We  give  some  examples  to  illustrate  some  important  points  with  respect  to  the  canonical  form. 

Example  I:  The  symbol  associated  with  the  Cauchy-Riemann  equations  is 


( 


*6 

*6 


(2.10) 


Its  determinant  — is  irreducible,  therefore,  this  matrix  is  not  C-reducible. 

The  factorization  of  the  determinant  into  polynomials  of  smaller  degree  does  not  guarantee  that 
the  corresponding  matrices  are  C-reducible. 


Example  II:  The  matrix 


3 


(2.11) 


whose  determinant 


*6 

*6 

*6 

0 

*6 

0 

(*6)(-^?  +  e2  +  ^3)  (212) 

is  a  product  of  two  lower  order  polynomials,  is  not  C-reducible.  This  can  be  seen  by  observing  that 
the  eigenvector  corresponding  to  the  eigenvalue  depends  on 


P-Reducibility 

For  numerical  applications  one  may  allow  a  more  general  transformation  in  reducing  a  system 
to  block  triangular  form  as  suggested  by  the  following  definition. 


Definition  2:  A  matrix  polynomial  P{C)  -'^oid  to  be  P-reducible  if  theiv  exist  invertible  matrix 
polyriomials  such  that  S(^)P(^)  =  A(^)T~’(^),  where  A(^)  is  a  matrix  polynomial  of 

block  tria7igular  form. 

The  matrix  given  in  Example  is  P-reducible.  This  can  be  seen  from  the  identity 


0 

0  \  1 

(  *6 

*6 

*6 

iiz 

*6 

0 

t-6  -*6  /  ' 

0 

Thus,  using 


*6 

1 

0  \ 

/  1 

0 

0  \ 

\e 

*6 

0 

0 

0 

0 

/ 

Vo 

*6 

/ 

(2.13) 


/  1  0  0  \ 

Si0  =  T-\0=  0  *6 

\  0  / 

we  get  a  reduction  to  a  block  upper  triangular  form  given  by 


(2.14) 


/  *6  1  0  \ 

A=  1^1'^  0  (2.1.5) 

\  0  0  / 

By  using  the  operators  corresponding  to  S,T  for  the  differential  problem  one  may  end  up  with 
other  variables  than  the  desired  ones.  This  reduction  is  still  important  for  analyzing  the  equation 
but  the  original  variables  may  need  to  be  solved  for  as  well.  One  can  define  an  extended  system  in 
which  the  canonical  form  is  combined  with  the  transformation  of  variable  and  get  a  block  triangular 
form  as 


j’-i 

0 


(2.16) 
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This  form  which  we  caU  the  extended  canonical  form  may  be  required  only  when  the  matrices 
S,T~^  depend  on  In  the  other  case  the  transformation  from  the  canonical  variables  to  the 
original  ones  is  done  by  an  algebraic  mapping.  In  the  applications  discussed  in  the  next  section  the 
extended  canonical  form  will  be  found  useful  for  three  dimensional  problem. 

Assuming  that  S,T  transform  the  matrix  polynomial  P{^)  into  block  triangular  form,  we  obtmn 
a  similar  reduction  for  the  differential  equation,  namely, 

{SP{D)T){T-^u)  =  Sf  (2.17) 

that  is,  in  terms  of  a  new  set  of  variables  T~^u  we  get 


AiD)v  =  Sf 

The  boundary  conditions  B{D)u  =  g  are  transformed  into 


(2.18) 


Bi,{D)v,  =  iQg)i  i=l,...,fc  (2.19) 

and  we  assume  that  for  each  j  —  I,. .  .,k  the  operator  Pjj{D)  with  the  boundary  condition  Bjj{D) 
form  a  well  posed  problem. 

The  simplest  set  of  boundary  conditions  are  those  for  which  the  boundary  operator  in  the 
transformed  variables  is  diagonal,  and  each  corresponding  subproblem  is  well  posed.  We  refer  to  a 
set  of  such  boundary  conditions  as  canonical. 

All  the  arguments  made  above  applies  to  linear  constant  coefficient  problems.  The  variable 
coefficient  case  and  the  nonlinear  case  may  have  triangular  forms  for  the  principal  part  only. 

It  should  be  mentioned  that  the  use  of  S,T~^  which  are  not  polynomials  can  reduce  F(^)  to 
upper  triangular  form.  However,  the  elements  of  this  form  are  not  polynomials,  therefore,  do  not 
correspond  to  differential  operators.  While  such  a  representation  is  still  useful  in  analysis  of  well 
posedness  of  the  boundary  value  problem,  it  is  not  as  useful  in  developing  numerical  methods  for 
the  problem.  The  representations  discussed  in  this  paper  involve  transformations  that  correspond 
to  differential  operators  and  can  be  used  in  numerical  implementation  of  certain  iterative  techniques 
that  are  based  on  these  forms.  These  issues  will  be  discussed  in  a  separate  paper. 


3  Canonical  Forms  for  Incompressible  Flow 

In  the  applications  we  present  in  this  and  the  next  section,  although  nonlinear,  a  full  reduction  to 
triangular  form  is  presented,  without  involving  any  linearization. 

The  steady  state  incompressible  inviscid  equations  in  nonconservative  formulation  [1]  are  given 

by 


divV  =  0 

(V  •  V)V  +  Vp  =  0 


(3.1) 


which  can  be  written  in  operator  form  as 


(3.2) 


/div  0\/v\/0\ 

y  Q  grad  )\p)  [  0  J 

where  Q  is  a  diagonal  matrix  operator  of  dimension  d  x  d  whose  diagonal  entries  are  the  scalar 
convection  operator  Q  =  V  -V. 

Analyzing  the  determinant  of  this  system  reveals  its  structure  which  is  necessary  for  analyzing 
possible  boundary  conditions  and  discretization  issues.  Denoting  the  matrix  operator  above  by 
L,nc,  and  freezing  the  coefficients  at  some  constant  flow,  Vq,  one  obtains 


cletL.„,(0  =  -l^r'{Vo-0  (3.3) 

where,  ^  =  (6i  •  • 

The  determinant  which  is  already  factored  into  irreducible  factors  reveals  important  properties 
of  the  system  of  the  incompressible  Euler  equations.  While  the  operator  A  is  elliptic,  the  operator 
Q  is  hyperbolic  with  respect  to  the  stream  direction.  This  system  is  therefore  of  a  mixed  type. 
Some  components  of  the  system  display  hyperbolic  behavior,  while  the  rest  show  elliptic  behavior. 
The  power  d  -  1  in  the  factorization  of  the  determinant  suggests  that  there  are  d  -  1  quantities 
that  are  governed  by  hyperbolic  subsystems.  This  forms  a  subsystem  of  the  full  equation  which  is 
weakly  coupled  to  the  rest. 

We  proceed  by  identifying  the  hyperbolic  components  of  this  system.  Using  the  relation 


(V  •  V)V  =  - V(V  •  V)  -  V  X  curlV  (3.4) 

and  defining 

P  =  P+\q^  =  (3.5) 

one  obtain  Crocco’s  form  of  the  equations  of  motion,  namely. 


div(V)  =  0 

-  V  X  curlV  +VP  =  0 


(3.6) 


Let  (61,62,63)  be  an  orthonormal  basis  such  that  61  is  in  the  direction  of  V,  and  (61,62,63) 
form  a  right  hand  system,  i.e.,  61  x  62  =  63,62  x  63  =  61,63  x  6]  =62 

Taking  the  inner  product  of  the  momentum  equations  with  V  and  using  V  •  (V  x  curlV)  = 
(V  X  V)  •  curlV  =  0,  we  obtain 


V- VF  =  0  (3.7) 

This  equation  implies  that  the  total  pressure  is  constant  along  streamlines.  Thus,  the  total 
pressure  P  is  one  of  the  canonical  variables  we  are  looking  for,  and  it  corresponds  to  a  factor  V  •  V 
in  the  factorization  of  the  system.  For  three  dimensional  flows  there  exists  another  quantity  that 
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is  governed  by  a  a  hyperbolic  operator.  From  physical  insight,  a  quantity  that  exi..ts  in  thiee- 
dimensiona]  and  not  in  tv^'o  dimensional  flows  may  have  to  do  with  vorticity.  The  vorticity 
defined  by 


oj  =  curlV  (3.8) 

will  be  shown  to  play  an  important  role  in  the  decomposition  of  the  system. 

It  can  be  easily  verified  that  the  vorticity  component  in  the  direction  of  the  velocity  cannot  be 
determined  from  the  momentum  equations.  However,  we  can  obtain  the  following 

—  Cj  •  (V  X  oj)  +  Cj  •  VP  =  0,  j  =  2,3  (3.9) 

which  gives 

-(ej  X  V)-w  +  ejVP  =  0,  j  =  2,3  (3.10) 

Using  the  definition  of  we  get  the  following  equations 

-7(62  •a>)  +  e3- VP  =  0 

7(e3-aj)  +  e2-VP  =  0  ^ 

which  can  be  interpreted  as  equations  that  determine  two  components  of  the  vorticity  vector. 

The  other  component  of  the  vorticity  is  obtained  from 


diva?  =  0  (3.12) 

which  is  a  compatibility  condition  (div  curl=0).  Decomposing  u>  as 


oj  =  fi\  +  w*’ 
w*-  •  Cl  =  0 


(3.13) 


and  substituting  (3.13)  in  (3.12)  and  using  the  continuity  equation  one  obtain  an  equation  for  (i, 
namely. 


V  •  V/3  +  divw*"  —  0  (3.14) 

Therefore,  I3  is  also  governed  by  a  hyperbolic  equation.  Note  that  while  P  is  constant  along 
streamlines,  (3  admits  a  more  complicated  structure  depending  on  the  behavior  of  P  in  cross-stream 
directions. 

Summarizing,  we  get  the  following 

divV  =  0 

curlV  —  u>  =  0 

V  •  V/3  -I-  div{u»-L)  =  0 

w-*- =  (62  •a>-‘-)e2 -t- (63  •u»-‘-)e3  (3.1.5) 

-7(62  •  oj-^-)  -1-  63  •  VP  =  0 
7(63  •  w-*-)  -I-  62  •  VP  =  0 
761  •  VP  =  0 
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which  can  be  written  in  a  matrix  form  as 


div 

\ 

—curl 

qe  1  e2  e3 

V  •  V  dive2  dive3 

0 

1 

*0 

0 

63  •  V 

0  0  q  ' 

62  •  V 

V  V  / 

[7) 


(3.16) 


Observe  that  although  the  upper  block  of  the  system  which  consists  of  of  four  equations 
(div,  cur/)  is  overdetermined,  its  solution  is  guaranteed  by  the  presence  of  the  fi  equation,  which 
states  that  divu»  =  0.  Another  way  of  looking  at  it  is  to  eliminate  the  divu>  —  0  equation  and  to 
regard  the  {div,curl)  system  as  an  equation  for  both  V  and  [i.  In  this  way  we  obtain  a  4  by  4 
system  for  four  unknown  quantities.  This  fixes  a  well  known  difficulty  of  solving  {dii\curl)  in  three 
space  dimensions. 

For  two-dimensional  flows  0  =  0  and  the  equations  take  a  simpler  form.  If  e),e2  lie  in  the  plane 
defined  by  the  flow  then  w*-  •  62  =  0.  Thus,  can  be  represented  by  a  scalar  quantity  which  we 
denote  by  w.  The  canonical  form  of  the  two  dimensional  incompressible  Euler  equations  reduces  to 


D.  Dy 

0 

0 

/  0  \ 

Dy  —  Dx 

-1  0 

V 

0 

0  0 

-9  ^Do 

U) 

0 

0 

0 

0  Q  y 

\^) 

where 

Do  =  vDj;  —  uDy  (3.18) 

These  representations  of  the  incompressible  inviscid  eruations  reveal  the  structure  of  solutions. 
That  is,  the  total  pressure  P  and  the  normalized  helicity  (with  respect  to  velocity)  evolve  along 
streamlines,  and  determine  the  vorticity  vector  field.  The  velocity  components  satisfy  an  elliptic 
system  of  a  (div,curl)  form. 

In  this  representation  which  we  call  the  canonical  form  of  the  incompressible  Euler  equations  the 
elliptic  part  has  been  separated  from  the  hyperboUc  part.  The  variables  {V,0,P)  will  be  referred 
to  as  the  canonical  variables. 

This  form  of  the  equations  suggest  a  natural  set  of  boundary  conditions,  which  we  refer  to  as 
canonical  boundary  conditions  for  the  corresponding  systems.  The  appearance  of  the  Laplacian 
operator  A  in  the  factorization  of  the  determinant  calls  for  one  boundary  condition  at  every  bound¬ 
ary  point.  The  term  implies  that  additional  d  —  1  boundary  conditions  are  to  be  prescribed 

at  inflow  points  of  the  boundary.  From  the  factorization  it  is  natural  to  prescribe  at  inflow  the 
canonical  variables  P,0,  and  the  extra  condition  to  be  imposed  at  every  point  is  V  •  n  where  n  is 
an  outward  normal.  Note  that  for  two-dimensional  problems  0  =  0  and  therefore  one  boundary 
condition  is  omitted. 
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4  Canonical  Forms  for  Compressible  Flows 

The  Euler  equations  in  conservation  form  [1]  in  terms  of  the  variables  {V,p,  H)  are 


div(pV)  =  0 

div(pV  M  V  +  p/)  =  0  (4.1) 

div{pV//)  =  0  =  0 

The  analysis  of  the  equations  is  simplified  if  one  moves  to  the  non-conservative  formulation, 

(  Qc  Qp  W  V  \  /  0  \ 

pQ  grad  0  j  p  =  0  (4.2) 

\0  0  pQ  )\h  )  \Q  ) 

where 

Qc  =  pdiv-H25^V.(V.V) 

Qv  =  (4.3) 

Qh  =  §tjQ 

Denoting  the  matrix  operator  above  by  Lg,  and  freezing  the  coefficients  at  (Vo,po, /fo)i  the 
determinant  of  the  symbol  gives 


det  Ui)  =  p^+2(Vo  •  0'((Vo  •  if  -  (4.4) 

where  cq  is  the  speed  of  sound  and  can  be  evaluated  using  the  relations  (where  subscripts  are 
omitted), 


s  =  p/p'T' 

H  =  +  W 

=  7p/p 


(4.5) 


The  factor  Vq  •  i  appearing  n  the  determinant  here  is  the  same  one  as  for  the  incompressible 
equation,  suggesting  that  there  are  d  quantities  that  propagate  along  streamlines.  The  rest  of  the 
system  is  governed  by  a  system  equivalent  to  the  well  known  full  potential  operator  whose 

character  depends  on  the  Mach  nui.  ber  M, 


=  q^/c^  (4.6) 

leading  to  an  elliptic  equation  for  M  <  1  and  hyperbolic  (in  the  stream  direction)  for  M  >  1.  As 
before  we  use  Crocco’s  form  of  the  equation  for  analysis  purposes. 


QoV  =  0 

—pV  XU}  —  pTVs  +  pV H  =  0 
div(pV//)  =  0. 


(4.7) 
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By  subtracting  H  times  the  continuity  equation  from  the  energy  equation  one  obtains 


pV  .  V//  =  0 

Taking  the  scalar  product  of  the  momentum  equation  with  V  leads  to 


(4.8) 


pV  ■  V// -pTV  •  V5  =  0  (4.9) 

Thus,  two  of  the  quantities  that  propagate  along  streamlines  are  the  total  enthalpy  H  and  the 
entropy  s.  As  the  determinant  suggest  there  is  one  more  quantity  in  three  dimensional  flow  that 
propagate  along  streamlines. 

To  find  that  quantity  we  follow  a  similar  path  as  in  the  incompressible  case  by  introducing  the 
vorticity  and  its  decomposition  as 


w  =  curlV 

uj  =  ^pV  +  u>-*-  (4.10) 

w-L  •  V  =  0 

The  equation  for  /?  is  obtained  the  same  way  as  in  the  incompressible  case,  i.e.,  using  divw  =  0, 
giving 


The  quantity  a?-*-  satisfies 


pV  •  Vfi  4-  divu)-^  =  0 


(4.11) 


-P9(e2  -UJ-*-)  -  pTea  •  Vs  +  pea  •  V/f  =  0  f4.]‘'>) 

pq{e3  •  uj-*-)  —  pTe2  ■  Vs  +  pe2  •  Vi7  =  0 

in  complete  analogy  with  the  incompressible  case. 

The  continuity  equation  is  simplified  by  subtracting  from  it  appropriate  multiples  of  the  energy 
and  the  entropy  equations,  given  QcV  —  0. 

Summarizing  we  get 


QcV  =  0 

curlK  —  <*>  =  0 

pV  •  V/J  +  div(w-^)  =  0 

uj-*-  -  (62  •u;-‘-)e2  -  (ea  •w‘-)ea  =  0 

-pq{e2  -w*-)  -  pTea  •  Vs  +  pea  •  V^  =  0 

pq(e3  -  w*-)  —  pT'e2-Vs  +  pe2  •  V/f  =  0 

-pTV  .  Vs  =  0 

pV  •  Vff  =  0 


or  in  matrix  form 
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/  Qc 

O 

II 

-curl 

pqe\  e2  e3 

pV  •  V  dive2  dive3 

0  -9  0 

0  0  p 

0 

—pTe^  •  V 
-pTei  ■  V 

0 

pes  •  V 
pe  2  •  V 

-pTV  •  V 

0 

v-v 

The  nature  of  the  subsystem 


(4.14) 


gcV  =  0  (4.15) 

curlV  =  a>  (4.16) 

can  be  studied  by  taking  the  gradient  of  the  first  equation  and  subtracting  from  it  the  curl  of  the 

second  equation.  Using  the  relation 

V^V  =  graddivV  —  curlcurlV  (4.17) 


and 

V  .  gradp  =  •  (V  .  VV)  (4.18) 

we  get 

V^V  +  grad(--7^^V  .  (V  •  VV))  =  -curfw  (4.19) 

pd(qy2) 

This  equation  is  elliptic  for  subsonic  flow  and  hyperbolic  with  respect  to  the  stream  direction  for 
supersonic  flows. 

In  two  dimensions  /3  =  0  and  the  canonical  form  of  the  Euler  equation  reduces  to 


where 


Di 

D2 

1  0 

0 

0 

-Dy 

D, 

-1 

0 

0 

0 

0 

-q 

0 

0 

0 

TpQ 

0 

0 

0 

0 

0 

pQ 

^  u  y 

f  0\ 

V 

0 

(jj 

0 

s 

\  H 

1  0/ 

D\  =  p/c^((c^  -  u^)Dx  -  uvDy) 
Di  =  /5/c^((c'^  -  v^)Dy  -  uvdJ) 
Dj  =  vDx  —  uDy 


(4.20) 


(4.21) 
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From  the  canonical  form  one  can  easily  see  that  a  natural  set  of  boundary  conditions  for  the 
steady  state  Euler  equations  is  the  following. 

At  subsonic  points  one  condition  have  to  be  specified  on  account  of  the  elliptic  part.  This  is 
the  compressible  analog  of  the  (div,curl)  system  for  the  incompressible  case.  This  condition  can  be 
pV  •  n  where  n  is  the  outward  normal  to  the  boundary.  A  compatibility  condition  has  to  hold,  that 
is,  the  integral  of  that  quantity  around  the  boundary  has  to  vanish,  similar  to  the  incompressible 
case.  At  inflow  boundary  points  d  more  conditions  has  to  be  specified  and  the  natural  choice  is 
as  these  quantities  are  carried  into  the  domain  by  the  convection  operator  Q.  At  outflow 
boundary  points  no  condition  is  required  for  the  operators  Q. 

At  supersonic  points  the  system  is  purely  hyperbolic  and  require  at  the  inflow  d  +  2  conditions, 
and  no  conditions  at  outflow.  In  that  case  it  is  natural  to  specify  the  quantities  (V,  H,s),  rather 
than  using  /J  for  one  of  them. 

5  Conclusion 

Canonical  forms  for  inviscid  flow  problems  have  been  derived.  In  these  forms  the  different  t>pes 
of  subsystems  have  been  separated  using  a  set  of  new  variables,  namely,  the  canonical  variables. 
These  forms  are  upper  triangular  operator  forms  in  which  the  building  blocks  of  the  systems  reside 
on  the  diagonal  and  the  interaction  of  the  different  type  of  subsystems  is  represented  by  the  off 
diagonal  blocks.  This  decomposition  of  systems  of  partial  differential  equations  allows  a  better 
insight  into  the  structure  of  solutions  which  is  essential  in  constructing  numerical  solutions. 

The  forms  described  here  suggest  new  iterative  solutions  of  the  inviscid  equations,  both  in  two 
and  three  dimensions.  These  will  be  discussed  elsewhere. 

The  different  approximations  used  in  fluid  dynamics  over  the  years  are  clearly  seen  in  the  canon¬ 
ical  form.  Assuming  the  flow  is  of  constant  total  pressure  P  one  obtain  a  (div,  cur/)  system,  which 
in  two  dimensions  is  nothing  but  the  well  known  Cauchy- Riemann  equations.  The  three  dimen¬ 
sional  implementation  of  a  (div,  cur/)  system  was  not  popular  since  the  overdetermined  system 
needs  to  be  discretized  carefully.  As  the  canonical  form  suggests,  the  introduction  of  0  into  that 
overdetermined  system  fixes  the  major  numerical  difficulties  related  to  existence  of  solutions. 

For  the  compressible  case,  by  assuming  that  the  total  enthalpy  H  and  entropy  s  are-constant, 
one  obtains  a  system  which  is  equivalent  to  the  full  potential  equation.  This  is  the  upper  left  block 
in  the  canonical  form. 
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